Collective vortex rectification effects in films with a periodic pinning array 
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The vortex ratchet effect has been studied in Al films patterned with square arrays of submicron 
antidots. We have investigated the transport properties of two sets of samples: (i) asymmetrical 
antidots where vortices are driven by an unbiased ac current, and (ii) symmetrical antidots where 
in addition to the ac drive a dc bias was used. For each sample, the rectified (dc) voltage is 
measured as a function of drive amplitude and frequency, magnetic field, and temperature. As 
unambiguously shown by our data, the voltage rectification in the asymmetric antidots is induced 
by the intrinsic asymmetry in the pinning potential created by the antidots, whereas the rectification 
in the symmetric antidots is induced by the dc bias. In addition, the experiments reveal interesting 
collective phenomena in the vortex ratchet effect. At fields below the first matching field (Hi), 
the dc voltage— ac drive characteristics present two rectification peaks, which is interpreted as an 
interplay between the one-dimensional motion of weakly pinned incommensurate vortex rows and 
the two-dimensional motion of the whole vortex lattice. We also discuss the different dynamical 
regimes controlling the motion of interstitial and trapped vortices at fields slightly higher than Hi 
and their implications for the vortex ratchet effect. 

PACS numbers: 74.78.Na., 74.40.+k, 05.40.-a, 85.25.-j 



I. INTRODUCTION 

Ratchets are systems with an asymmetric periodic 
potential capable of promoting unidirectional transport 
of classical or quantum objects subject to an unbi- 
ased fluctuating drive. For that purpose, the exter- 
nal forcing must be a non-equilibrium one, like an al- 
ternating force or a colored noise, since useful work 
can not be extracted from equilibrium fluctuations only. 
The ratchet effect was first used to explain biochem- 
ical mechanisms such as the intracellular transport of 
macromoleculesii Since then, parallel to advances in the 
theoretical understanding^ several devices have been 
proposed and fabricated in order to test the ratchet idea 
at the classical and quantum levels and to use the ratchet 
effect for particle segregation and motion rectification of 
electrons and fluxons^ Most of this work, however, is fo- 
cused on ratchet systems consisting of single objects or 
an assembly of weakly interacting objects. Much less is 
known about the collective rectification of strongly in- 
teracting systems via the ratchet effect and work on this 
subject carried out so far is mostly theoretical. 

Vortices in superconductors with a periodic array of 
pinning centers have proven to be an ideal system in 
which to study collective dynamical phenomena of par- 
ticles moving in a periodic potential. Several interesting 
dynamical effects have already been studied in these sys- 
tems. For instance, as demonstrated theoretically^ and 
experimentally£*£*£, periodic pinning potentials are able 
to guide vortex motion into high symmetry directions of 
the pinning lattice. Another remarkable dynamical ef- 



fect in these systems is the presence of Shapiro steps in 
the voltage-current characteristics when an rf current is 
coupled to the dc drive £ The steps appear due to phase- 
locking between the rf signal and the vortex motion over 
the periodic pining structure. 

Controllable motion of vortices by means of the ratchet 
effect was recently demonstrated in experiments per- 
formed on superconducting samples with periodic arrays 
of asymmetric pinning centersSalfi and for asymmetric 
configurations of symmetrical pinning sitesi These struc- 
tures are able to break the symmetry of the vortex pin- 
ning potential and promote vortex motion rectification 
when an unbiased ac current is applied. In Ref. |3, for 
instance, the vortex ratchet effect was demonstrated in 
a Nb film where the asymmetrical pinning potential was 
given by an array of triangular magnetic dots. In Ref.llCt 
we observed pronounced rectification effects in Al films 
where the symmetry of the pinning potential was bro- 
ken by engineering a composite configuration of antidot 
lattices with big and small antidots placed close to each 
other. 

Despite a great deal of theoretical effortsii*i2ii&±2ii£ 

aimed at the description of the ratchet dynamics of vor- 
tices in asymmetrical pinning potentials, these systems 
are far from being completely understood. The vortex 
motion itself, in the presence of strong pinning sites, is 
still a subject of intense study. Differently from many 
other systems, vortices are soft objects; their core is able 
to deform to better adjust to the available pinning po- 
tential. Due to this property, the usual approximation of 
treating vortices as localized overdamped particles may 
break down and the resulting vortex ratchet effect can 
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be quite different from, for instance, the effects known 
for overdamped Brownian ratchets. Indeed, as we have 
recently shown, the dynamics ruling the voltage rectifica- 
tion in a superconducting film with an array of asymmet- 
rical antidots can be described by underdamped equa- 
tions of motion^ This means that this system can be 
regarded as an example of inertia ratchet systems, even 
though the true mass of a vortex is negligible. 

In the present paper, we investigate in details vortex 
motion rectification in Al films with nanoengineered an- 
tidot arrays as a function of ac excitation, magnetic field 
and temperature. Two systems are investigated: (a) 
square arrays of asymmetrical antidots excited by a sym- 
metrical ac current, an intrinsic ratchet system, and (b) a 
square array of symmetrical antidots excited by an ac cur- 
rent coupled to a dc bias, i.e. a tilted-potential ratchet. 
For fields below the first matching field, Hi = ^o/flp 
(where $o — h/2e is the flux quantum and a p is the an- 
tidot lattice constant), the vortex response in both sys- 
tems is characterized by hysteresis, which suggests under- 
damped motion, and a double rectification peak at some 
incommensurate fields. These results are interpreted in 
terms of a minimal inertia ratchet model and molecular 
dynamics simulations of "inertial" vortices. The simu- 
lation results indicate that the double peak observed at 
incommensurate fields can be explained in terms of plas- 
tic deformation of the vortex lattice. For fields higher 
than Hi, our data reveal that the motion of interstitial 
vortices is essentially overdamped, in contrast with the 
unusual underdamped dynamics of the vortices trapped 
by the antidots. This result suggests that the underlying 
inertial effect is connected to the strong interaction be- 
tween the vortices and the antidots and not to the real 
mass of the vortices. 

The paper is organized as follows. In Sec. [H] we de- 
scribe the sample preparation, the experimental proce- 
dure and the experimental results. In Sec. lIIII we present 
the details of our underdamped model and the molecular 
dynamics simulations. The experimental data, as well 
as a comparison between modelling and experiment, are 
discussed in Sec. IIVI Finally, in Sec. [V] we draw our 
conclusions. 



II. EXPERIMENTAL RESULTS 

A. Sample preparation and experimental details 

The experiments where performed on three different 
Al thin films, thermally evaporated on a SiC>2 substrate 
patterned by electron beam lithography. For all samples 
the pattern was designed in a cross shaped geometry to 
allow four point measurements in the two perpendicu- 
lar directions [see Fig. El (a)]. The cross consists of two 
300-/um-wide strips containing the nanoengineered array 
(period a p = 1.5 /im) of pinning sites. This gives a value 
of 0.92 mT for the first matching field. The first two 
samples (AAD1 and AAD2) have the same unit cell con- 



TABLE I: Superconducting parameters for the three samples 
studied. 





AAD1 


AAD2 


SAD 


Thickness t (nm) 


38 


37 


22 


Critical temperature T c (K) 


1.469 


1.438 


1.532 


Mean free path l e (nm) 


3.0 


8.1 


3.4 


Coherence length £(0) (nm) 


59 


97 


63 


Penetration depth A(0) (nm) 


225 


136 


210 


GL parameter k 


3.83 


1.40 


3.32 



sisting of a small (300 x 300 nm 2 ) and big (600 x 600 
nm 2 ) antidot separated by a thin superconducting wall 
(90 nm) and approximate thicknesses (38 and 37 nm). 
The third sample (SAD) has a symmetric unit cell con- 
sisting of 600 x 600 nm 2 square antidots and a thickness 
of 22 nm. Figs. ^ (b) and (c) show atomic force micro- 
graphs (AFM) of samples AAD2 and SAD, respectively. 
The dashed lines depict the unit cell (1.5 x 1.5 /im 2 ) of 
the respective sample. The white dots observed in the 
images are particles of resist layer which did not come off 
during the lift off procedure. These particles are innocu- 
ous to our experiment. 

Table 0] shows the superconducting parameters for all 
three samples. The superconducting critical tempera- 
ture T c was obtained by using a resistance criterion of 
10% of the normal state resistance in a zero-field mea- 
surement. The resistive transition for all samples was 
sharp, with a width typically of AT C ~ 3mK . From the 
residual resistivity at 4.2 K, the electron mean free path 
l e was found. From this we have calculated the coher- 
ence length £(0) and the penetration depth A(0) using 
the dirty limit expressions. All samples are type-II su- 
perconductors, i.e., their Ginzburg-Landau parameter is 
k = A(0)/£(0) > 1/V2. Note that, although samples 
AAD1 and AAD2 share the same pattern and have simi- 
lar thicknesses, their mean free paths were quite different, 
which is a result of the different evaporating conditions 
(base pressure and growth rate) in which the samples 
were prepared. 

The transport measurements were performed in a cryo- 
stat using ac and/or dc currents applied at the end of 




(a) (b) (c) 



FIG. 1: Layout of the Al samples, (a) Cross shaped geometry 
of the samples. AFM images of sample AAD2 (b) and SAD 
(c). The dashed lines indicate the corresponding unit cell. 
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the cross legs in the cc-direction. The magnetic field was 
applied perpendicularly to the xy plane. In all ac mea- 
surements, the applied ac current was a sinusoidal wave, 
i.e., lit) — I ac sin(27r/t), where / is the excitation fre- 
quency and I ac is its amplitude. The voltage was ac- 
quired in four point mode by an oscilloscope, used as an 
analogue-digital converter. The waveform was averaged 
over several cycles (typically 200-500 cycles) in order to 
enhance the signal-to-noise ratio. The dc voltage was 
then obtained by integrating the "clean" waveform. 



B. Rectification by unbiased asymmetrical antidots 

Fig. 0(a) shows a contour plot of the dc voltage across 
the sample {Vdc) as a function of reduced magnetic field 
(H/H\) and current amplitude (I a c) at a temperature 
T = 0.98T C and excitation frequency / = 1 kHz. At 
H = 0, one can observe a clear sign inversion of the dc 
voltage. Note that the voltage drop due to vortex motion 
is given by Vdc = L-((v)xB), where B is the flux density, 
(v) is the time averaged vortex velocity, and L is the 
distance between the voltage contacts. Thus, the change 
of Vdc sign as H (and, consequently, B) changes its sign 
means that vortices and antivortices are rectified to the 
same direction. This simple assumption demonstrates 
unambiguously the intrinsic rectifying properties of the 
asymmetrical antidots configuration. A similar behavior 
of Vdc{H, I ac ) was also observed in sample AAD1 and its 
plot for the positive field side was presented in Ref. 0. 
Another important effect demonstrated by Fig. [21 (a) is a 
clear enhancement of the rectified voltage near the first 
matching field (see also Ref. Fig. 2). This can be 
explained by the high symmetry of the vortex lattice at 
this field range which induces coherent vortex motion, 
that is, all vortices contribute to the rectified voltage in 
a constructive way. 

In panels (b) and (c) of Fig. [2 we show detailed plots of 
the Vdc- Iac characteristics of samples AAD2 and AAD1, 
respectively. Similarly to single-object ratchet systems, 
these curves are characterized by a rectification window, 
in which the rectified voltage increases monotonically, 
and by the presence of a tail, in which Vdc drops smoothly 
towards zero. The rectification window is defined by two 
critical forces, Fi, which marks the onset of rectification, 
and F2, for which rectification is maximum. These forces 
also determine the effective asymmetry of the pinning 
potential, which we define here as a = 1 — F1/F2. For 
some field values (usually, for fields close to Hi or Hi/2) 
a second rectification peak is observed, suggesting an in- 
terplay between two ratchets (c.f. Sec. I1V() . each having 
their own critical forces [see Fig. Hfc)]. However, a clear 
second peak was not identified in the measurements on 
sample AAD2, which is probably due to a lack of res- 
olution in our measurements in the narrow rectification 
window of this sample. 

For fields higher than H\ , a reversal in the net vortex 
motion is observed at lower amplitudes [see the Vd c (I a c) 
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FIG. 2: (a) Contour plot of the dc voltage measured across 
sample A AD 2 as a function of reduced field and ac ampli- 
tude at T — 0.98T C , and frequency of 1 kHz. (b) Detail of 
the amplitude dependence of the dc voltage shown in (a) for 
fields H/Hi = 0.95, 0.97 and 1.03. (c) Ac-drive amplitude 
dependence of the dc voltage measured across sample AAD1 
at T = 0.976T C and / = 1 kHz for fields H/Hi = 0.49, 0.82 
and 0.96. 

curve at H/H 1 = 1.03 in Fig.^b)]. In Ref.H a similar 
voltage reversal was observed at fields above 3Hi , which 
was the saturation field of the sample studied. The au- 
thors interpreted this result by arguing that the vortices 
trapped by the triangular magnetic dots arrange them- 
selves in such a way that the interstitial vortices expe- 
rience an effective asymmetric potential with opposite 
asymmetry. In our experiment, the saturation field was 
Hi. This can be inferred from the strong drop of the 
critical current just above this field, which indicates that 
weakly pinned interstitial vortices come into play. Then, 
in fact, the voltage reversal observed at fields higher than 
Hi can also be interpreted as a reversal in the recti- 
fication of interstitial vortices. However, in our case, 
the reverse asymmetry is not induced by an asymmet- 
ric arrangement of vortices inside the pins, but rather by 
the current distribution (encircling the asymmetric an- 
tidot configuration) created by a single-quantum vortex 
trapped in the double-antidot pins. 

By comparing the results obtained on samples AAD1 
and AAD2, one can notice the striking difference in the 
rectified voltage of both samples. For the similar reduced 
temperatures, the maximum dc voltage of AAD2 was sev- 
eral times smaller than that of AAD1. This can be ex- 
plained by the fact that the coherence length of AAD2 
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0.25 




ac Amplitude (mA) 

FIG. 3: Dc voltage measured across sample AAD2 versus 
ac amplitude for H/Hi = 0.95 and / = 1 kHz at different 
temperatures, as indicated. Inset: Plot of the maximum dc 
voltage at this field as a function of temperature. 



at temperatures close to T c may be bigger than the anti- 
dots (for instance £,aad2 (T = 0.97T C ) ~ 560 nm). In this 
situation, the vortices are not able to resolve the details 
of the antidots and the effective asymmetry is strongly 
suppressed. This is very clear by calculating the asymme- 
try ratio of both samples at a relatively low temperature 
T = 0.976T C (where thermal escape can be neglected) 
and reduced field H/H\ — 0.96, which gives a w 0.09 for 
AAD1 and a w 0.02 for AAD2. 

The coherence length, can also play an important role 
in the temperature dependence of the rectified voltage. 
Fig. shows Vdc-Iac characteristics of sample AAD2 at 
several temperatures for the same reduced field H/H\ = 
0.98. The maximum rectified voltage seems to increase 
monotonically until a temperature T ~ 0.982T C , above 
which it starts to drop, resulting in a bell-shaped T de- 
pendence (see inset). In Brownian ratchets, a decrease 
in the dc response as the temperature increases is in- 
deed expected»i& The reason is that, when the Brownian 
motion of particles is enhanced by thermal noise, more 
energy is needed to move them in a given direction. How- 
ever, in the vortex case near the critical temperature, 
increasing temperature induces not only vortex fluctua- 
tions but also a strong increase of the vortex size, which 
leads to a decrease in the effectiveness of the pinning 
potential and its asymmetry. At low temperatures, the 
reduction of fluctuations and vortex size should lead to a 
monotonic increase in rectification. Nevertheless, when 
vortices are small enough and their fluctuations are weak, 
they are more sensitive to the background pinning caused 
by the natural inhomogeneities of the sample. Due to the 
disorder induced by the background pinning, the anti- 
dot pinning potential may become less effective at lower 
temperatures, thus causing a reduction of the rectified 
voltage. 

In Fig.0Ja)-(d) we show the voltage output waveform 
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FIG. 4: (a)-(d) Time evolution of the voltage output of sam- 
ple AAD1 at a temperature T = 0.973T C , field H = 0.96-Hi, 
and ac amplitudes I ac = 477 (a), 502 (b), 544 (c) and 704 /xA 
(d). (e),(f) Voltage output of sample SAD at T = 0.966T C , 
H — 0.82J/i, a dc bias — 64 fiA, and ac amplitudes: I ac = 
553 and 569 /iA (e); I ac = 624 and 680 fiA (f). (g),(h) Volt- 
age output of sample SAD at T = 0.966T C , H = 1.50f/i, a dc 
bias Idc = 64 fiA, and ac amplitudes: I ac — 102, 221 and 292 
fiA (g); and I ac = 320, 413, and 479 fiA (h). 



in one cycle of the ac drive for different drive amplitudes 
of the Vdc-Iac characteristic presented in Fig. |2f c) (open 
circles). Note that, as discussed above, here negative 
voltage corresponds to positive vortex velocity, since the 
applied field is positive. Panels (a) and (b) correspond 
to points in the rectification window, where the ratchet 
system behaves as a half-wave rectifier, whereas (c) and 
(d) correspond to points in the rectification tail. At each 
half-period of the ac drive, the corresponding half-period 
of the voltage output in all cases (a-d) is asymmetric; vor- 
tices seem to be depinned at a given force and repinncd 
at a lower force value, giving rise to hysteretic behavior. 
Such a hysteresis was also observed in sample AAD2, but 
somewhat smoothed due to the strong effects of fluctua- 
tions in this sample. As we shall see below, this indicates 
that the system is ruled by underdamped dynamics. 



C. Rectification by dc biased symmetrical antidots 

Since the symmetrical antidot configuration of sam- 
ple SAD does not produce rectified vortex motion by it- 
self, we externally induce an asymmetry in the system 
by "tilting" the pinning potential with a small dc cur- 
rent bias. In Fig. [5Ja), we show the contour plot of the 
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FIG. 5: (a) Contour plot of the dc voltage measured across 
sample SAD as a function of reduced field and ac amplitude 
at T = 0.966T C , frequency of 1 kHz, and a dc bias of 68 
fiA. (b) Detail of the amplitude dependence of the dc voltage 
shown in (a) for fields H/Hi = 0.43, 0.82, and 1.09. The 
dashed line indicates the saturation voltage V sa .t = R n Idc = 
2.8 mV. (c) Ac-drive amplitude dependence of the dc voltage 
at T = 0.973T C , / = 1 kHz, H/Hi = 0.93, and different tilts, 
Ida =130, 46 and 22 fiA. 



dc voltage across sample SAD in the H-I ac plane at a 
temperature T = 0.966T C for a dc bias of Idc = 64 /iA 
and ac excitations of 1 kHz. [Note that the steps observed 
in this plot are just artifacts induced by the low resolu- 
tion of the field sweep (0.05 mT).] In this experiment, no 
flip of the dc voltage sign is observed. This is indeed ex- 
pected because here the polarity of a vortex does matter 
in determining its net velocity. For instance, for a posi- 
tive current bias, i.e., in the positive x direction, a vortex 
with positive polarity experiences an positive slope in the 
pinning potential (due to the Lorentz force pointing to 
the negative y direction), whereas an antivortex senses a 
negative slope. Thus, as the magnetic field changes its 
sign, so does the asymmetry of the pinning potential, in 
such a way that the cross product (v) x B keeps its sign. 

In Fig. 02b), we present a detailed plot of the Vdc-Iac 
characteristics at fields H/Hi = 0.43, 0.82, and 1.09, 
for the experiment of Fig. [SJa). Here, in comparison 
to the Vdc-Iac characteristics of the asymmetric antidot 
samples (Fig.|2Jl, the important difference is that the rec- 
tification tail approaches asymptotically a finite voltage 
value, which depends on the applied dc tilt, whereas the 



rectification tail of the intrinsically asymmetric samples 
reaches zero voltage at high ac drives. In the rectification 
window, however, samples SAD and AAD1 behave in a 
similar way. For instance, at fields close to Hi, the sec- 
ond rectification peak is also seen in the Vdc-Iac curves 
of sample SAD. 

Fig. Efc) presents the rectification effect in sample 
SAD for different dc bias, I dc = 22, 46, and 130 ^A at 
H/Hi = 0.93 and T/T c = 0.973. At this field, a double 
rectification peak is observed for all dc biases. The dis- 
tance between peaks seems to be constant, whereas the 
rectification window is wider for higher tilt values. The 
window is roughly defined by the critical values I c — Idc 
and I c + Idc, where I c = 501 /iA is the critical depinning 
current obtained from conventional dc transport mea- 
surements. Thus, the system is given an asymmetry of 
a — 2Idc/{I c + Idc)- For the tilt Idc = 22 fj,A, this gives 
a rs 0.086, which is similar to the intrinsic asymmetry 
observed in sample AAD1 at T = 0.973T C . Another note- 
worthy point is that the main rectification peak decreases 
relatively to the saturation voltage for higher tilts. This 
is consistent with the fact that, in the limit of the critical 
tilt {Idc = I c ), no rectification peak should be observed 
because there the energy barrier for the easy flow direc- 
tion disappears completely. In this limiting case, Vd c {Iac) 
is expected to grow monotonically from zero to the sat- 
uration voltage. 

The time evolution of the voltage output for differ- 
ent drive amplitudes of the Vdc-Iac characteristics shown 
in Fig. |Sfb) is presented in Fig. Qfe)-(h). We consider 
two different fields, H/H x = 0.82 (e and f) and 1.50 (g 
and h), corresponding, respectively, to filled circles and 
crossed circles of Fig. [Sfb) . In the H/H\ = 0.82 case, the 
asymmetrical shape of the half-periods of V(t) resembles 
that observed in sample AAD1 [shown in panels (a)-(d) 
of this figure], that is, here a hysteresis in the depinning- 
repinning process is also observed. The main difference 
is that in the high drive limit [panel (f)], the amplitude 
of the positive half-cycle is higher than the amplitude 
of the negative one, due to the extrinsic nature of the 
asymmetry in sample SAD, whereas for sample AAD1 
the amplitude of both half-cycles are approximately the 
same. 

For fields higher than Hi [Fig. 0] (g) and (h)] , there ex- 
ists an amplitude range for which the voltage half-cycles 
are symmetrical, that is, no hysteresis is observed. Since 
this behavior is only present at low ac amplitudes, we 
can assume that the symmetrical waveforms correspond 
to motion of interstitial vortices. At HjH\ = 1.5, the 
response is strongly enhanced in the positive half-cycle 
above an amplitude I ac = 292 /iA, suggesting that at this 
value the vortices trapped by the antidots start moving. 
The voltage waveform corresponding to this amplitude 
and above are asymmetrical [panel (h)], resembling the 
waveforms observed at fields below Hi. A similar evo- 
lution of the voltage waveforms as the ac amplitude is 
varied was also observed in the samples with asymmetri- 
cal pinning sites for fields H > Hi. 
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III. THEORETICAL MODELLING 
A. Inertia ratchet model 

Our experimental data clearly demonstrate a pro- 
nounced hysteresis in the voltage as a transport current is 
cycled. Typical systems exhibiting pronounced hystere- 
sis in their transport properties are driven underdamped 
particles in periodic potentials. In these systems, the 
hysteresis in the depinning-repinning process is a conse- 
quence of the finite inertial mass of the particles, which 
has important effects on their transport properties. For 
instance, in microscopic boundary lubrication problems, 
the depinning of the (inertial) lubricant particles from the 
periodic potential of a solid substrate is usually identified 
as the static friction force whereas the repinning of the 
particles is identified as the dynamical friction forced 
Other "inertia" ratchet systems are asymmetrical under- 
damped Josephson junction arrays^ and dc SQUIDs 19 . 
In these two examples, the rectified object is the phase 
of the superconducting order parameter across a junction 
and the inertia is given by the shunt capacitance of the 
junctions. This problem is formally the same as that of 
the mechanical inertia ratchet. 

The rectification of inertial particles in a ratchet 
potential has been studied theoretically by several 
groups pI&iiLSiiSi The dynamics of inertia ratchets is more 
complex then their overdamped counterpart due to the 
existence of both regular (for moderate underdamped) 
and chaotic (for extremely underdamped) dynamics. In 
the regular regime, the amplitude dependent dc response 
presents a strong enhancement in the rectification win- 
dow and a faster decaying tail, as compared to the char- 
acteristics of overdamped particles. Inspired by the fact 
that, in our experiment, (i) the rectified voltage demon- 
strates hysteretic behavior and (ii) the Vd c -I a c character- 
istic generally has pronounced rectification effect in the 
rectification window, we propose a description of vor- 
tex dynamics in antidot systems in terms of a simplified 
underdamped model, in which vortices (still treated as 
particles) are given an apparent mass to account for the 
observed hysteresis. 

We first consider the simpler case of a single vortex 
in a ratchet potential, here representing the asymmetric 
antidot array. The equation of motion is given by 

Mx = -r)x - Up(x) + Asin{vt), (1) 

where 77 is the viscous drag coefficient, x is the particle 
position. The pinning potential U p (x) is modelled in a 
simplified way as the double-sine function: 

/ s r [2ttx\ (3 . /4vra;\l 

u >w = - ^HvH sm Ur)J (2) 

The parameter [3 defines the asymmetry of the pinning 
force (a = 2/3/(1 + /3), for (3 < 1/4), the critical forces 
being given by F w = 2ir(l — f3)U/a p , for drive along the 
(weak) positive x direction, and F s = 2ir(l + (3)U/a p , 
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FIG. 6: Time averaged particle velocity as a function of 
ac drive amplitude calculated in the single-particle inertia 
ratchet model [Eq. for the viscosity-to-mass ratios rj/M = 
0.1, 2.5, 4 and 10 (in units of f p ). Inset: Contour plot of aver- 
age velocity as a function of amplitude and rj/M; the dashed 
lines define the rectification window in this case (a = 0.095). 



for drive along the (strong) negative x direction. Note 
that, since here we do not take into account thermal 
fluctuations, the critical forces which define the rectifi- 
cation window of this system are given deterministically 
by F\ = F w and F 2 = F s . 

The main characteristic frequencies setting the limits 
of the dynamical regimes of Eq. Q are the libration fre- 
quency in the potential minima, which for a < 0.2 is 

given approximately by f p = yJu/Ma^, and the inverse 

relaxation time 7 = rj/M. The quasi-adiabatic regime 
(i.e., close to the zero- frequency adiabatic limit), corre- 
sponds to frequencies much smaller than these two char- 
acteristic frequencies. In the numerical solution of Eq. 
we have chosen excitation frequencies / < 10 -4 / p and 
7 > 0.1/p, corresponding to the quasi-adiabatic regime. 

In Fig. we show the net dc velocity versus ac drive 
amplitude ((v)-F ac characteristics for different values of 
7. Forces are in units of 2irU/a p and frequencies (and 
7) in units of f p . For these runs, we have chosen the 
asymmetry factor a — 0.095, which is close to the typical 
values of J1//2 observed in sample AAD1. The inset 
shows a detailed diagram of (v) as a function of A and M. 
Notice the crossover between the spread rectification in 
the extreme overdamped limit and the sharp rectification 
in the extreme underdamped regime. For rj/M > 4/ p , the 
curves are essentially overdamped, following the behavior 
of the adiabatic solution of a deterministic overdamped 
ratchet^ For 77/M < 2f p , the curves are much more 
pronounced in the rectification window and have a faster 
decaying tail. We shall compare these results with the 
experimental data in Sec. IIVI 
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B. Molecular dynamics simulations 

Now, we consider the 2D problem of a lattice of "iner- 
tial" vortices interacting with a 2D array of asymmetric 
pinning sites, which is modelled as two interpenetrating 
square sublattices of strong and weak attractive poten- 
tials. The sublattices have spacing a p (—1 here) and are 
shifted with respect to each other by a distance d = 0.2a p . 
All pinning centers are modelled by Gaussian potential 
wells with a decay length R p . Hence, the pinning force 
exerted by the k-ih pin of one sublattice on the z-th vor- 
tex reads 



T71S.UJ 



-F, 



R„ 



■ exp 



where represents the position of the z-th vortex and 
R s k ' w is the location of the k-th pin of the weak (w) and 
the strong (s) sublattice, respectively. The intensity of 
the individual pinning force is denoted by F SiW . We de- 
fine the ratio between the weak and strong pinning forces 
as e = F w /F s . In our simulation, forces (per unit length) 
are taken in units of F = 4>§/87r 2 A 3 , with A the ef- 
fective superconducting penetration depth. Vortices are 
considered to interact logarithmically, the corresponding 
repulsive vortex- vortex force being given by 



^vv ) F v 



„o ^2 Arj 



where = (r^ — rj)/\vi —Tj\. F vv q denotes the strength 
of the vortex- vortex interaction. The underdamped equa- 
tion of motion for vortex i is given by 



Mr, = F L {t) + ~F vv (r^ ) + P p (rO - m 



(4) 



The ac driving Lorentz force is F^(t) = Asin(ujt), where 
oj = 2ir/P and P is the period of the ac drive. r\ is the 
viscosity coefficient and M is the effective mass of an 
individual vortex. This equation is integrated by molec- 
ular dynamics (MD) simulations. Parameters used in the 
simulations are R p = 0.13ao, A = 0.26ao. F vv q — 0.06, 
F p s = 0.5 and e = 0.925, time step t = 1/150. All 
the frequencies of the ac sine- wave drive are smaller than 
10~ 6 /to, i.e., into the quasi-adiabatic regime. For the 
mass and viscosity we have chosen, respectively, M = 1 

and r) = 1, which, in terms of f p — «/ U p /Ma p (here 

Up = F s R p /2 is the amplitude of the individual strong 
pinning potential), leads to a viscosity-to-mass ratio 
rj/M = 2.44/p. The MD simulation results will be dis- 
cussed in the next section. 



IV. DISCUSSION 

In Fig. [7| we present a fitting of our inertia ratchet 
model to the experimental data of samples AAD1 and 



SAD. The data is plotted in normalized units. The nor- 
malized critical forces used in the model were the same 
observed in the respective measurement. In this way, 
the only adjustable parameter is the viscosity-to-mass 
ratio, 7 — rj/M, which is given in units of the libra- 
tion frequency f p . For comparison, we have shown in 
Fig. the Vdc-Iac calculated close to the overdamped 
limit 7 = 10. As it is clearly seen from this graph, the 
usual overdamped model fails to describe the ratchet ef- 
fect of vortices in an antidot lattice, more specifically in 
our Al samples AAD1 and SAD. On the other hand, the 
underdamped model with 7 values close to the under- 
damped limit provides very good fitting to our experi- 
mental data. For the measurement of sample AAD1 at 
T = 0.973T C and H/H t = 0.99, for example, we used 
a = 0.095 (which gives approximately the same normal- 
ized critical forces observed in this experiment) and the 
best fit was obtained with 7 = 1.8 [see Fig. [7| (a)-inset], 
which is close to the underdamped limit. For the sample 
SAD at H/Hx = 0.98 and T = 0.966T C , we used a = 
(i.e. a symmetrical sinusoidal potential) and added to 
the driving force a dc component F^ c = Q.0QQF p (where 
F p = 2ix\J I a v is the critical force of the pinning poten- 
tial), in agreement with the Idc/Ip ratio observed in the 
experiment. In this case the best fit was obtained with 
7 = 2.1 [Fig. 0(b)]. 

The success of an underdamped model to describe the 
voltage induced by vortex motion is not a trivial re- 
sult. Vortex dynamics is usually treated in terms of over- 
damped particle models because it is a well established 
fact that in most samples the vortex mass is negligi- 
ble compared to the damping coefficient ^ Nevertheless, 
these models do not consider the effect of vortex deforma- 
tion that can lead to an irreversibility in the depinning- 
repinning process and, consequently, to an apparent in- 
ertial effect. In recent Ginzburg-Landau simulations of 
a superconductor with an antidot lattice, it was shown 
that the process of depinning of a vortex from an antidot 
can be characterized by strong current-induced elonga- 
tion of the vortex core cross-section, which results in a 
hysteresis in the pinning mechanism^ This effect seems 
to mimic the behavior of an underdamped particle in a 
periodic potential and can be responsible for the hystere- 
sis observed in our experiments. 



A. Vortex rectification for H < H\ 

In our experiment, the Vd c -I ac curves for some fields 
below Hi are characterized by a double rectification peak. 
As we pointed out in our previous publication, the sec- 
ond rectification peak can be explained in terms of plas- 
tic deformation of the vortex lattice siS In other words, 
the existence of a plastic dynamical phase with a lower 
critical current gives rise to the appearance of a second 
ratchet system, which competes with the main ratchet, 
given by the motion of the whole vortex lattice. The 
consequences of the existence of two ratchets with dif- 



8 



o 
> 
o 

"D 



1.0 


" (a) 




1.0 

! 0.8 




|H = 0.99H 1 | 


0.8 








Y c 




0.6 






it 0.2 
'^0.0 


j 




0.4 










0.2 




ii 

i i 


III IV 

, !i , i 




1.2 1.4 


0.0 


Li' 


|h = 

I 7— 


0.96 



±2 u L i J i I i I . I ■ I i I 

1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 







0.6 
0.4 
0.2 
0.0 



- 1 I 1 I 1 I 1 I 


1 1 1 1 1 






Experiment 






Num. Fit 


















H = 0.98H 1 | 




1 . 1 1 . 1 


i.i. 





0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 

Normalized ac Amplitude 

FIG. 7: (a) Main frame: normalized dc voltage versus cur- 
rent amplitude (circles) of sample AAD1 at T = 0.973T C 
and H = 0.96-Hi and numerical fit obtained from the dou- 
ble ratchet model described in Sec HVf see also Sec. with 
7 = 1.9 and \ — 0.5. The dashed vertical lines define different 
dynamical regimes of the system (see text). Inset: same as 
the main frame but with H — 0.99-Hi. The fitting (line) was 
obtained integrating Eq. Qwith 7 = 1.8. (b) Normalized dc 
voltage versus current amplitude (circles) of sample SAD at 
T = 0.966T C and H — 0.98-f/i using an excitation frequency 
of 1 kHz and a dc bias of 46 fj,A (which gives a reduced tilt. 
The fitting (line) was obtained using 7 = 2.1. 



amplitude F ac for reduced fields HjH x = 1, 0.94, 0.92, 
and 0.88. The simulations reproduce quite well the dou- 
ble rectification peak observed experimentally. Panels 
(e) and (f) show the time dependence of the center-of- 
mass (cm.) velocity of the vortex lattice at vortex den- 
sities H/Hi = 1 and 0.92, respectively. As in our ex- 
periments, here the waveforms of vortex motion during 
each half-period are asymmetrical, i.e., present hystereti- 
cal behavior, which is a result of the finite apparent mass 
of the vortices. In particular, the simulated waveforms 
at H/Hi = 0.92 [panel (f)] evolve similarly to the ex- 
perimental data of sample AAD1 presented in Fig.|3|(a)- 
(d). In both, experiment and simulations, four different 
ratchet phases (I-IV) are clearly identified [these phases 
are defined by vertical dashed lines in the V& c -l a c charac- 
teristic of Fig. [7] (a) and in the (v)-F ac characteristics of 
Fig. 0(b) and (c)]. In phase I, only a fraction of vortices 
move in the easy direction of the pinning potential. No 
motion to the hard direction is observed. In phase II, the 
whole vortex lattice is half- wave rectified, i.e., all vortices 
move to the easy direction and no vortices move to the 
hard one. In phase III, we have motion of all vortices in 
the easy direction and motion of a fraction of them to 
the hard direction. The averaged rectified voltage drops 
abruptly at the onset of motion to the hard direction but 
increases again at slightly higher amplitudes, since the 
main ratchet (the whole vortex lattice) is still in its rec- 
tification window, giving rise to the second rectification 
peak. Finally, in the phase IV, all vortices move back 
and forth, but in an asymmetric fashion, resulting in the 
long tail observed in the experiments and simulations. 



ferent critical forces can be hcuristically illustrated by 
our one particle underdamped model. Let us consider 
two independent inertia ratchets. The first one repre- 
sents the whole vortex lattice, with critical forces F w and 
F s . The second ratchet represents the fraction x of the 
vortices that deform plastically and has critical forces 
F w — F w an( i F' s < F s . In this way, the total dynamics 
is simply given by the solution v(t) of the first ratchet 
for F(t) > F w and F(t) < —F s , and, otherwise, by the 
solution v'(t) of the second ratchet, normalized by the 
factor x- This simple model provides an excellent fitting 
to the experimental data for H < H\. For instance, to fit 
the data of the experiment performed in sample AAD1 
at H = 0.96^1 and T/T c = 0.973, we solved Eq. Q for 
both ratchets independently, using for the critical forces 
the normalized values obtained experimentally, and ad- 
justed the parameters 7 and x- In this case, we obtained 
7 = 1.9 and x = 0-5 [see Fig. 0(a)]. 

To get a better insight on how such plastic dynamics 
can lead to a double peak in the Vd c -I a c characteristic, 
we performed MD simulations of "inertia!" vortices for 
fields H < H\. The parameters used in the simulations 
are given in Sec. IIIIH 3. Fig. |S] (a)-(d), show the simu- 
lated net velocity of vortices (v) as a function of drive 



To illustrate the interplay between the plastic dynam- 
ics and motion of the whole vortex lattice, we plotted 
in Fig. [S] (g)-(i) the vortex trajectories for the simula- 
tion at H/Hi = 0.94 and a few drive amplitudes. At 
F ac — 0.2 (g), all vortices are still pinned (just oscillating 
inside the pinning centers). Note that, since the number 
of vortices does not match the number of pinning sites, 
there is a distribution of vacant pins, that is to say, of 
discommensurations in the vortex lattice. In this case, 
and also for H/Hi = 0.92, there are some vortex rows, 
with a denser distribution of discommensurations, which 
depin more easily. Thus, at F ac — 0.205 (h), only the 
weakly pinned vortex rows are dragged away, giving rise 
to a one-dimensional channelling dynamics characterized 
by a strong shear of the vortex lattice. Finally, at a 
higher drive, F ac = 0.210 (i), the potential barriers for 
the strongly pinned vortex rows are also overcome and all 
vortices move to the easy direction of the pinning poten- 
tial. At lower fields, such as H/H\ = 0.88 shown in panel 
(d) of this figure, the distribution of vacancies is approx- 
imately homogeneous in such a way that all vortex rows 
has approximately the same depinning force. Thus, no 
channelling dynamics is observed and the (v)-F ac char- 
acteristic at this field presents a single rectification peak. 
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FIG. 8: Molecular dynamics simulations 
results for 77 = 1 and M = 1. Pan- 
els (a)-(d): net dc velocity as a function 
of the ac driving force amplitude in nor- 
malized units (see text) at reduced fields 
H/Hi = 1 (a), 0.94 (b), 0.92 (c), and 0.88 
(d). The dashed lines in (c) and (d) define 
different dynamical regimes of the system 
(see text) . Panels (e) and (f ) : time evolu- 
tion of the cm. velocity at H/Hi = 1 (e) 
and 0.92 (f). Panels (g)-(i): vortex tra- 
jectories (lines) for H/Hi — 0.94 at drive 
amplitudes 0.2 (g), where vortices are still 
pinned, 0.205 (h), where part of the vor- 
tices move as ID channels, and 0.21 (i), 
where all vortices move. Vortices are rep- 
resented by the small dots and the double 
pinning sites are represented by the dou- 
ble discs. 



B. Vortex rectification for H > Hi 

At fields higher than Hi , the interpretation of vortex 
motion in terms of underdamped dynamics may break 
down because of the presence of interstitial vortices. 
Since they do not interact directly with the antidots, their 
dynamics is expected to be similar to that of vortices in a 
plain film, that is, essentially overdamped. This was in- 
deed suggested by the amplitude dependence of the volt- 
age waveforms V(t) [see Fig.|2|(g) and (h)]: whereas for 
high amplitudes the signal in a half-cycle is asymmetric, 
at lower amplitudes, the signal may be perfectly sym- 
metric, that is, completely reversible. Since apparently 
hysteresis only come out when the vortices trapped in the 
antidots are set into motion, one may then conclude that 
the apparent inertia comes from the interaction between 
vortices and antidots rather than from the real vortex 
mass. 



V. CONCLUSION 

In conclusion, we have studied the vortex ratchet effect 
in superconducting films with nanoengineered arrays of 
asymmetric and symmetric antidots. In the asymmetric 
antidot sample, we observed rectification of an unbiased 



ac current due to the intrinsic asymmetry of the system, 
thus demonstrating that this is an intrinsic vortex ratchet 
system. In the symmetric antidot sample, the rectifica- 
tion effect was due to the dc tilt of the periodic (symmet- 
ric) potential induced by the applied dc bias. The vor- 
tex rectification is strongly dependent on temperature 
and magnetic field. The temperature dependence indi- 
cates that the coherence length plays an important role 
in determining the effective asymmetry produced by the 
asymmetric antidot array; when the vortex size becomes 
comparable to the big antidot size, the rectification is 
strongly suppressed. 

The field dependence of the rectification effects found 
in our experiments provides the following observations in- 
dicating collective behavior in the ratchet dynamics: (i) 
rectification is enhanced near the first matching field, ev- 
idencing the importance of vortex-vortex interactions in 
the ratchet effect (at the first matching field, the vortex- 
vortex interactions cancel out by symmetry) ; (ii) at fields 
above Hi , we observed a sign reversal in the ratchet effect 
at low amplitudes, which indicates that weakly pinned 
interstitial vortices are rectified to the opposite direction 
due to the inverse asymmetry produced by the current 
distribution of the vortices trapped in the antidots; and 
(iii) at some fields below Hi , a double rectification peak 
was observed, which can be explained by the competition 
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between rectification of (weakly pinned) incommensurate 
vortex rows and the rectification of the whole vortex lat- 
tice. 

In addition, our data revealed that the ratchet effect 
for fields lower than H± is dominated by underdamped 
dynamics. These results were corroborated by molecular 
dynamics simulations of an underdamped vortex ratchet 
model. We argued that the inertial effect actually comes 
from vortex deformation due to the strong interaction 
with antidots, rather than from a finite vortex mass. In- 
deed, at fields higher than H±, we observe that the dy- 
namics of interstitial vortices is essentially overdamped, 
since these vortices are not pinned directly by the anti- 



dots. 
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